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1 Introduction 



Hadronic atoms represent new tools for probing strong interaction dynamics at low en- 
ergies [1, 2, 3, 4, 5, 6]. When the hadronic constituents are pseudoscalar mesons, then 
the probe concerns the properties of spontaneous chiral symmetry breaking in QCD. The 
simplest representative in that case is the inr atom (pionium), which was produced and 
studied in the DIRAC experiment at CERN [7] and to which many theoretical studies 
were devoted [8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. At the next level in that 
category one finds the irK atom, the properties of which are tightly related to those of the 
SU(3) x SU(3) chiral symmetry breaking. In particular, the energies and widths of the 
atomic levels depend on the strong interaction nK scattering lengths and through them 
on various order parameters of chiral symmetry breaking. The theoretical interest of those 
quantities justifies the preparation of new experimental projects for the production and 
study of irK atoms [21]. On theoretical grounds, the properties of the ttK atom were 
recently studied in detail by Schweizer [22] using the approach of nonrelativistic effective 
field theories to the bound state problem [23]. 

The present work is devoted to the study of the properties of the nK atom using the 
quasipotential-constraint theory approach [24, 25, 26, 27, 28] which was also applied to 
the pionium case [15] and which consists of a relativistic and covariant three-dimensional 
formulation of the bound state problem. Our aim is to calculate the 0(a) corrections to 
the lowest-order formulas for the energy shift and decay width of the 7rK atom, where 
a is the fine structure constant, and to provide the means of extracting with sufficient 
precision the values of the strong interaction S'-wave irK scattering lengths from future 
experimental data. 

We introduce with respect to our previous approach to the pionium problem a slight 
modification in that we formulate from the start the bound state problem almost on the 
mass shell. This is suggested by the substantial simplification one gains in order to reach 
the final results. The leading nonrelativistic formulas, which are of order a 3 , are expressed 
mainly in terms of threshold properties of the strong interaction on-mass shell scattering 
amplitudes. In general, bound state problems in quantum field theories are formulated 
in terms of off-mass shell scattering amplitudes or kernels; on the other hand, effective 
lagrangians, like the one used in chiral perturbation theory, give rise to a proliferation 
of terms, some of which do not contribute on the mass shell. Therefore, one expects 
many cancellations among unphysical quantities, which are automatically realized in a 
formalism based from the start on the use of on-mass shell scattering amplitudes and 
their minimal analytic continuations. 

However, because of the presence of infrared divergences, the on-mass shell formalism 
ceases to be consistent at higher orders in QED. Generally, up to 0(a> 4 ) effects in the bound 
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state problem infrared divergences are unambiguously recognized, isolated and subtracted 
or cancelled, according to the method of approach. This is an implicit consequence of the 
Ward identities (and their analogues for mass terms) satisfied by the theory. At higher 
orders, for instance in the treatment of the Lamb shift problem, infrared divergences can 
no longer be isolated without ambiguity in an on-mass shell formalism of the bound state 
problem. The recourse to an off-mass shell formalism becomes compulsory at this stage 
[29, 30]. For the present problem, since the precision that is sought does not necessitate 
the evaluation of 0(a 5 lna _1 ) effects, this difficulty will not show up. 

The main corrections, of order a, to the lowest-order formulas can be represented by 
three groups of terms: i) the pure elecromagnetic corrections, arising beyond the Coulomb 
potential; ii) electromagnetic radiative corrections to the strong interaction scattering 
amplitudes, including isospin symmetry breaking effects; iii) second-order perturbation 
theory effects in the perturbative expansion of the bound state energy. The sum of those 
corrections is found to be of the order of a few per cent relative to the leading terms. Our 
results agree with those obtained by Schweizer [22]. 

The plan of the paper is the following. In Sec. 2, we present the general bound state 
formalism that we use. In Sec. 3, we consider the specific case of the irK atom, for 
which we adopt a coupled channel formalism. In Sec. 4, we calculate the lowest-order 
contributions to the real and imaginary parts of the energy shift. In Sec. 5, pure QED type 
corrections are evaluated. In Sec. 6, electromagnetic radiative corrections to the strong 
interaction scattering amplitudes are taken into account. In Sec. 7, a general property 
of cancellation of divergences of three-dimensional integrals, present in our formalism, is 
displayed. Sec. 8 deals with second-order effects of the perturbation theory expansion 
of the bound state energy. A summary of results follows in Sec 9. Some evaluations of 
integrals are presented in the Appendix. 

2 Bound state formalism 

It is generally recognized, on the basis of hamiltonian formalism, that relative times 
and relative energies of particles of multiparticle systems should not play a dynamical 
role in relativistic theories. Constraint theory [31, 32] allows, through the use of first- 
class constraints, the elimination of redundant variables, respecting at the same time the 
symmetries of the theory (in the present case the Poincare invariance). For a two-particle 
system, described by momenta p\ and P2, with physical masses mi and 77*2, the following 
constraint 

C(P,p) = (jp\ - (ml -mj) = 0, P = Pl+P2 , p=-(p 1 -p 2 ), (2.1) 
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eliminates the relative energy in a covariant way. That constraint also respects the sym- 
metry between the two particles and remains valid on the mass shell or in the free case. 

We next consider a prototype quantum field theory with a given lagrangian describing, 
among others, the interaction between two spin-0 particles, 1 and 2. Let T be the on- 
mass shell elastic scattering amplitude of the process 1 + 2 — > V + 2'. [T is defined from 
the S matrix with the relation Sfi = 5fi + (27r) 4 <5 4 (P — P')Tfi, where P and P' are the 
total momenta of the ingoing and outgoing particles.] T is then a function of the two 
Mandelstam variables s = P 2 and t = q 2 , with q — (jp\ — p^). For convenience, we define 
a modified scattering amplitude T: 

f(s,t) = ^=T(s,t). (2.2) 

Since the equations we are developing are manifestly covariant, the total momentum 
P of the system defines a priviledged direction and all momenta can be decomposed 
along longitudinal (one-dimensional) and transverse (three-dimensional) components with 
respect to P. However, to simplify notation, we shall henceforth stick to the center-of- 
mass frame (P = 0) and represent vectors with their temporal and spatial components. 

The starting point of the present formalism is the postulate that T satisfies, by means 
of an effective propagator g Q , a three-dimensional Lippmann-Schwinger type equation 
leading to the definition of a kernel or a potential V(s,t): 

V = f-Vg f. (2.3) 

This equation can be used to calculate V either iteratively from T, or from a perturbative 
expansion of T itself. The iterative integration is three-dimensional taking into account 
constraint (2.1). Thus, if the two particle momenta linking T to V are (pi—k) and (p2 + k), 
total momentum being conserved, constraint (2.1) applied to them yields k = and 
integration concerns the three-momentum k. The latter integration does not, however, 
preserve the mass-shell character of T; it forces the corresponding momentum transfer 
— k 2 to take unphysical values down to — oo. The situation here is very analogous to that 
of nonrelativistic dynamics. We consider it as corresponding to a minimal extension of 
the scattering amplitude to the off-shell, in the sense that it does not imply introduction 
of new terms, but simply extension of the domain of validity of the expression of T to a 
larger one. 

The expression of g is chosen so that V is hermitian in the elastic unitarity region [26] . 
Noticing that constraint (2.1) implies equality of the Klein-Gordon operators of particles 
1 and 2, 

H (s, p) = {p\ - m\)\ c = {pi - m 2 2 )\ c = b 2 (s) - p 2 , 
„ . . s 1 . „ „. fm? — ml') 2 



g is chosen as the propagator associated with these operators: 

, , 1 1 , . 

9o{s,p) = — —, r = 7777 7, — . (2.5) 

On the mass shell, one has p 2 = &q( s ) and Hq vanishes. However, inside the integration 
domain of the iterative series, p is replaced by (p — k) and H takes the value (2p.k — k 2 ). 

It is advantageous in some instances to rewrite the scattering amplitude in terms of a 
two-particle irreducible kernel. Introducing the two-particle irreducible kernel K and the 
two-particle free propagator Go, one has T = K + KGqT and T takes the form: 

f = K + KG f, (2.6) 

where K = iK/(2\fs) and Go = —i2y/sGo, with 

Go = GioG2o, G a o = — — — — , a = 1,2. (2-7) 

Pa m l + 

Then the potential V can be expressed in terms of the kernel K: 

V = K(l-(G -g )Ky\ (2.8) 

the integrations involving Go being four-dimensional. 

We assume that V has been calculated (exactly or approximately) from T or K. 
Because of the mass shell condition imposed on the external particles, it is, like T, a 
function of s and t only: V = V(s,t). In x-space, obtained by Fourier transformation 
with respect to q (= (p — p')), it is a local function of x; all the nonlocal character of 
the interaction is now contained in the energy dependence of V. That feature does not 
remain true in an off-mass shell formulation, where V may depend on the operator p 2 . 

We next introduce the Green function associated with V and T (integrations on in- 
ternal variables are implicit): 

G( S ;p,pO = (2ii) 3 5 3 (p-p')g (s,p) + g (s,p)V(s,-k 2 )G(s;p-k,p') 

= 9 + 9 Tg . (2.9) 

Since G is an off-mass shell quantity, the external momenta p and p' are not restricted 
in the above equations to the mass shell and T is extended off the mass shell through the 
continuation of the variable t. 

In the presence of bound states, G develops poles in s. (We assume C, P and T 
invariances.) Let Sq be the mass squared of such a bound state (s > 0); then, in the 
vicinity of sq, G behaves as: 

G~ — , (2.10) 



S^SQ 



(s - s + is) 
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where \P' is the adjoint of Since (by construction) the kernel V does not have singu- 
larities, at least in the vicinity of the two-body threshold, one deduces from Eqs. (2.9) 
and (2.10) the wave equation 



So" 1 " V 



# = o. (2.ir 



The normalization and orthogonality properties of wave functions are obtained in the 
standard way [25, 27, 33]. Let us assume that the spectrum contains several nondegenerate 
bound states, labelled by an integer n, with masses squared s n > 0. (Generalization to the 
degenerate case is straightforward.) Designating by G~ l the inverse of G [G^ 1 = g$ l — V] 
and writing the wave equation of the bound state n as G~ 1 {s n )'^ n = 0, one easily arrives 
at the equation 

\ s s n J s s n 

Taking successively the limits s — > s n and s — > s m (m ^ n), one obtains the normalization 
and orthogonality conditions: 

r d*p A( dG-\sh _ , s 

r d 3 p T+ (G-\s n ) -G- l (s m )\ T 

The perturbative calculation of energy shifts from a zeroth-order approximation is also 
obtained with standard techniques [27, 33, 19]. Let us assume that the kernel V can be 
separated into two parts, 

V = V 1 + V 2 , (2.15) 

so that the solutions corresponding to V\ are known. Let G\ be the Green function 
associated with V\\ 

G 1 = g + g ViG 1 . (2.16) 

We designate by <p n the corresponding bound state wave functions with masses squared 
sffl . The complete Green function is constructed from G\ and the kernel V-i- 

G = G X + G 1 V 2 G. (2.17) 

The scattering amplitude due to V2 is obtained from V2 and G\. 

f 2 = V 2 + V 2 G 1 f 2 , (2.18) 

from which one deduces: 

G = G 1 + G 1 T 2 G 1 . (2.19) 
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To obtain the perturbative expansion around the bound state n, one isolates the corre- 
sponding pole in G\. 

Gi = G[ + ^ . (2.20) 
{S — Sn + ie) 

The reduced scattering amplitude T' 2 is constructed from the reduced Green function G[ 
and V 2 : 

f' 2 = V 2 + V 2 G' l f! 1 . (2.21) 
The relationship between T 2 and T' 2 is: 

f 2 =(l-f' ) % (2.22) 

V (s — sh J +ie)J 

Expressing in Eq. (2.19) G\ and T 2 in terms of G[, T' 2 and the pole term, one obtains the 
decomposition (integrations are implicit) 

G = G [ + G[ % G[ + (1 + G^) + m) , (2,3) 

from which one deduces the value of the mass squared of the bound state n [19]: 

s n = s^ + (^ n f^ n ). (2.24) 

Expansion of T 2 ' in terms of V 2 and G[ according to Eq. (2.21) yields the perturbative 
series of s n . Up to second order in V 2 , the expression of s n is [27, 33]: 



4 0) + \ U v 2 vO + (rf v 2 g; v 2 + v 2 vofat ^) 



(2.25) 

JO) v 7 



3 7ri^ system 

The 7T.K" atom is formed in the charged sector (ir~K + or its charge conjugate) under the 
effect of the Coulomb interaction and decays under the effect of the strong interaction 
predominantly into the neutral sector (tt°K or its charge conjugate). Branching ratios 
of other decays, involving photons, do not exceed a fraction of a per cent. It is natural 
to treat the ttK system by means of a coupled-channel formalism by generalizing the 
formalism developed in Sec. 2 with a matrix notation. 

We label with the index c the quantities related to the charged sector (n~K + ) and 
with the index n those related to the neutral sector (n°K ). Because of the decay process 
tt~ K + — > 7T°K°, the energy of the bound state becomes complex with a negative imaginary 
part. The scattering amplitudes and Green functions involving the above sectors have a 
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common pole at the position of the complex energy of the bound state. We introduce a 
two-component wave function \1/ as: 




*=UJ (3 ' 1) 

and define the potential V in matrix form in the corresponding space: 

V = ( Vcc Vcn | . (3.2) 




The iteration effective propagator g [Eqs. (2.3), (2.4) and (2.5)] is now composed of 
two propagators: 

, (!; ), *- = (f °X (3.3) 

g 0c and g Qn are defined with the physical masses of the particles, with respective energy 
factors bl c (s) and bl n (s) [Eq. (2.4)]. 

The wave equation (2.11) takes now the form of two coupled equations: 

(9oc-Vcc)^c-V c A = 0, (3.4) 
-KA+te-Kn)^ = 0. (3.5) 



The wave function represents an outgoing wave created by the charged state; it can 
be eliminated in favor of \l/ c , yielding the wave equation for the latter wave function: 

g£V c = V cc $> c + V cn (l - g 0n Vnn) _1 3onKc*c. (3.6) 

This is the bound state equation describing the properties of the irK atom. 

The potentials V are calculated from the Lippmann-Schwinger type equation (2.3), 
written now in matrix form in terms of the scattering amplitudes of the processes tt~K + — > 
ir~K + , ir~K + — > 7T°/i , 7r°/i° — > 7T~ K + , 7i°K° — > 7T° K° , which we designate respectively 
by T cc , T nc , T cn and T nn , T being defined in Eq. (2.2). The relationships of the components 
of V and T are: 

^cc = T cc Vcc 9oc Tcc V C n 9on Tnc i 
V C n T cn V C c 9oc Ten V C n 9 On Tnn i 

V nc T nc V nc gQ C T cc V nn gQ n T nc , 

V n n T nn V nc go c T cn V nn gQ n T nn . (3.7) 

When electromagnetism and isospin symmetry breaking are switched-off, one remains 
with the strong interaction or hadronic amplitudes T h in the isospin symmetry limit; 
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these are related to the isospin invariant amplitudes T 1 in the following way: 



T — T - (t 1 / 2 — T 3 / 2 

^ cn,h -L nc,h o V 



f nn , h = \{f^ + 2f^). (3. 



3 

One also defines the isospin even (+) and odd (— ) amplitudes: 

T+ = I(f 1 / 2 + 2T 3 / 2 ), 

T- = i(T 1 / 2 -T 3 / 2 ). (3.9) 
The decomposition into partial waves is done according to the formula (in the cm. frame) 

oo 

-2v/iT 7 (s, t) = 16vr ^(2^ + l)4(s)P^(cos 0), (3.10) 

1=0 

while the scattering lengths and effective ranges are defined from the threshold expansion 
of the real part of t\: 

TZetj(s) = ^(p 2 r(4 + ^p 2 + 0((p 2 ) 2 ). (3.11) 

The iteration series (3.7) defining the potentials involve with the presence of the effec- 
tive propagators g three-dimensional diagrams, which we call constraint diagrams due to 
the constraint (2.1) that is used there (see Sec. 2). They cancel the s-channel singularities 
of the scattering amplitudes in the scattering region and provide potentials that are real 
and regular in the total energy variable, being thus appropriate for a continuation to the 
bound state region. 

For Feynman diagrams involving QED parts, the calculation of V can be done with 
a simultaneous perturbative expansion of the scattering amplitude in the fine structure 
constant a. Since the bound state region is close to the two-particle threshold, one can use 
for the evaluation of the magnitude of the corresponding terms the threshold expansion 
method developed by Beneke and Smirnov [34] , which consists of making the expansion 
already at the level of the integrand, using dimensional regularization, by recognizing the 
various types of infrared singularities that might arise. The momenta that are relevant 
for that type of analysis are classified as potential, soft, ultrasoft and hard [34]. 

In order to apply perturbation theory to Eq. (3.6) for the evaluation of the bound 
state energy levels, it is natural to choose the nonrelativistic Coulomb potential as the 
zeroth-order potential, associated with the nonrelativistic kinetic energy. To this end, 
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we separate from the total energy Pq the mass term by defining (in the cm. frame) the 
binding energy S: 

^r s = P = m ^_ +mK+ + £, ( 3 . 12 ) 

Expanding s with respect to £ , one has for the cm. momentum squared factor &o c (s) [Eq. 
(2.4)]: 

,2 / X _ o..C fl , l m V + m K+ - m n- m K+\c\ .. m *- m K+ 



b' 0c (s)~2pL£ 1 + ^ £ , /i=- 7r ~ v (3.13) 

The Coulomb potential, which we designate by V^, appears in the nonrelativistic expan- 
sion of the potential V cc with the coefficient 2[i: 

n 

V cc = 2 t xV c + V cc , V c = --. (3.14) 

Therefore, we can divide the whole wave equation by 2/i and recover at zeroth order the 
nonrelativistic hamiltonian. The quadratic part in 8 in the expression of b^(s) above can 
then be treated as a part of the perturbation. Furthermore, by absorbing in the definition 
of wave functions the factor (jj/(m w - + m^+)) 1//2 (within the present approximation) one 
recovers from Eq. (2.13) the usual nonrelativistic normalization for the zeroth-order wave 
functions <p n : 

1 ti<Pn = I- (3.15) 



(2tt)3™ 

Thus, if satisfies the Schrodinger equation 

E^^={^ + V c )^ £i 0) = -f^ n>£+l, (3.16) 

with n, i and m representing the principal, orbital and azimuthal quantum numbers, 
respectively. 

Similarly, isolating, as in Eq. (3.12), from the bound state energies the mass term, 
by writing ^fs^ = m w - + m K+ + £ n and \J sffl = m^- + m K + + we recover from 

Eq. (2.25) the perturbation series of the nonrelativistic theory, written here up to second 
order: 



£n = ^i 0) + (^){(riV 2 ^) + (riV 3 GiV aVn 



ml- + m 2 K+ - m 7T -m K+ 
2m 7r _m^+(m 7r _ + m K+ ) 



S=E£ 

E^ )2 . (3.17) 



The last term comes from the relativistic corrections of the left-hand side of the wave 
equation (3.6) (cf. Eq. (3.13)). The last term of Eq. (2.25) has been discarded, since it 
does not contribute to order a 4 . 
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The bound state energy S n is in general complex. We designate by E n its real part; 
its imaginary part is equal to — T n /2, where Y n is the decay width of the bound state into 
two neutral mesons. We thus have: 

i 

£ n = E n - -Y n . (3.18) 

The Green function G\ [Eq. (2.16)] that is associated with the zeroth-order potential 
Vi, is now essentially the Green function of the Coulomb potential, which we designate 
by G c , the expression of which has been given by Schwinger [35]: 

(2tt) 3 5 3 ( P - p') 1 4vra 1 



G c (E,p,p') = 2fiG, 



(E - pV(2/i)) (E - pV(2/,)) (p - p') 2 (E - p'V(2/i)) 

1 l 
■ 4Trar]I(E, p, p') 



(£- P 7(2/i)) ' v (S-p«/(2/*)) 
- 2/x|y 0c + g Qc 2pV c g Qc - g 0c 8ir[iar]Ig 0c y (3.19) 

where / and rj are defined as 



I(E, p, p') = f dpp-v Up - p') 2 + (1 - p) 2 ^(E - p 2 /(2f,))(E - p l2 /(2pi)) 1 ^ 
Jo 1 cr 



a 



V = 1 = ■ (3-20) 
-E/(2fi) 



The quantity G' c will then correspond to G c from which the pole term around which 
perturbative expansion is organized is removed: G' c = 2fiG[ [Eq. (2.20)]. 

The potential Vi that appears in the perturbative formula (3.17) is obtained from the 
wave equation (3.6), by subtracting from the kernel of the right-hand side the Coulomb 
part [Eq. (3.14)]: 

V 2 = V cc + V cn (l - g, n V nn y l g, n Vn C . (3.21) 

In order to use this expression in Eq. (3.17), we first expand the factor (1 — g^Vnu)" 1 
in g 0n and evaluate the order of magnitude of each term of the expansion. The integration 



of g Qn yields the factor — iyoo n (s)/(47r) [Eq. (A. 5)], where the total energy is fixed at the 
bound state energy and b^ n {s) is defined in Eq. (2.4) with masses of the neutral mesons, 
m n o and m K o; it is equal to the square of the cm. momentum, which we represent by 
p*, of the system ir°K after the decay of the bound state: bl n (s) = p* 2 . The latter 
is essentially determined by the mass differences between charged and neutral mesons. 
Defining Am as 

2Am = m K+ - m K o + - m n o = 0.622 MeV, (3.22) 
one has approximately 

p* \ 2 Am 

-0.0029, (3.23) 
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which is a quantity of the order of a/2. Therefore, each g 0n can be estimated numerically 
as a quantity of the order of a 1//2 . Furthermore, the real part of the energy of the bound 
state receives contributions from even powers of g 0n , while the imaginary part (the decay 
width) receives contributions from odd powers of g Qn . With these estimates, Eq. (3.17) 
becomes: 



Here, AS n represents the energy shift with respect to the nonrelativistic Coulomb bound 
state energy (3.16). We have retained, for the real part of the energy, terms contributing 
to the orders a 3 and a 4 , while for the imaginary part of the energy (the decay width), we 



order terms are numerically negligible. 

4 Energy shift and decay width in lowest order 

We evaluate in this section the lowest-order contributions, 0(a 3 ) and 0(a 3 (Am//i) 1 / 2 ), 
to the real and imaginary parts of the energy, respectively; these come from the first 
two terms of the decomposition (3.24). Pure electromagnetic interactions beyond the 
Coulomb potential contribute only at 0(a A ) and hence can be ignored at the present 
level. A similar conclusion holds also for the interference terms between strong and 
electromagnetic interactions. Therefore, one has to consider in the potentials V cc , V cn 
and V nc solely contributions coming from strong interactions. 

All qualitative properties that are derived subsequently depend only on the short-range 
nature of the strong interaction, or equivalently, on the absence of massless particles in 
it; the particular model used for representing the strong interaction scattering amplitude 
is not relevant, although we will particularly refer to the chiral effective lagrangian [36, 
37, 38] as a prototype theory which will also be used for the numerical calculations. 

It is preferable here to evaluate V from its relationship with the two-meson irreducible 
kernel K, Eq. (2.8), which must be considered in its matrix form in the two-channel 
space. The kernel K is made of vertices and eventually of loops and may have momentum 
dependences. By analyzing, with the threshold expansion method [34], its behavior when 
integrated in Eq. (2.8) with g Q near the two-particle threshold (three-dimensional inte- 
gration), one finds that the eventual three-momentum dependences of K produce after 




have retained terms contributing to the orders o^Am/^) 1 ' 2 and a 4 (Am//i) 1 ^ 2 . Higher 
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integration momenta in the form p 2 or q 2 multiplying momentum independent pieces. 
Such terms have additional o? contributions with respect to the momentum independent 
pieces when considered in the QED bound states and thus can be neglected. Concerning 
the four- dimensional integration with Go, one notices that the presence of loop momenta 
in the numerator improves in general the infrared behavior of the integral; in this case 
it is only the hard momenta, which feel the ultraviolet behavior of the integral, that are 
concerned by their presence. Therefore, without loss of generality, the infrared behavior 
of integrals involving the kernel K can be studied by considering the latter as a constant 
and its value fixed at the two-meson threshold. In princple, the total energy is fixed at the 
bound state energy, but since K is a smooth function of the energy and since the binding 
energy is of order a 2 , one can continue without harm the energy of K up to threshold. By 
convention, one chooses for the strong interaction masses in the isospin symmetry limit 
the charged meson masses. 

With the above simplification, the integrations in Eq. (2.8) involve only Go and go- 
Here, however, the effective propagators g 0c and g Qn are defined with the meson physical 
masses [Eq. (2.1)]. It is then necessary, to carry out a consistent calculation, to also 
consider the propagators in Go with the corresponding physical masses. The details of 
the integrations are presented in the Appendix [Eqs. (A. 5), (A. 6) and (A. 7)]. A typical 
one-loop diagram, involving Go, and its constraint diagram, involving g , is presented in 
Fig. 1. 




Figure 1: Two- meson one- loop diagram and its constraint diagram, denoted with a cross. 



The integration of g 0c , at the bound state energy, provides the factor y — 6q c (s)/(47t) 
[Eq. (A. 5)]; it cancels a similar factor present in the integral of Go c [Eqs. (A. 6)- 
(A.7)]. Similarly, the inegration of g 0n , at the bound state energy, provides the factor 
~ ^\Aon( s )/(47r); it cancels the imaginary part of the integral of Go n . The remaining 
parts of the integrals of Go c and Go n are real smooth functions of the energy and can 
be replaced, dropping 0(a 2 ) corrections, by their values at the charged meson threshold. 
The above calculations can be repeated with all loops present in Eq. (2.8). One thus 
recovers the real part of the scattering amplitude at the charged meson threshold. Once 
the contributions of the effective propagators g have been taken into account, and in 
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order to avoid double-countings, one must take again in the remaining parts of the strong 
interaction amplitudes the equality of the masses of isospin partners, fixed at the charged 
meson masses. The effect of mass differences is taken into account separately within the 
calculations of electromagnetic and isospin symmetry breaking contributions. 

In summary, the strong interaction potentials that contribute at first-order pertur- 
bation theory to the complex shift of the bound state energy are given by the (charged 
meson) threshold values of the real parts of the strong interaction scattering amplitudes: 



V cc,h 
Vcn,h 



Vcc,h : 
Vnc,h 



Tie T CCjh 
■ He T cn h 



thr. 



thr. 



Tie T nc h 



thr. 



(4-1) 



where the subscript h refers to the purely hadronic part of the corresponding quantity. 
The threshold values of the real parts of the scattering amplitudes are given by the 5"-wave 
scattering lengths [Eqs. (3.8), (3.10) and (3.11)]. One finds: 



1/2 



3/2x 



47T 

Vcc,h = — — (2a 

47T i- 

Vcn,h = Vnc,h = ~^-y^[O'0 



1/2 



,3/2>, 



(4.2) 



Since those potentials are constant in momentum space, they yield delta-functions in 
x-space and their diagonal matrix elements become proportional to the square of the 
modulus of the wave function at the origin. In that case, only S-wave states contribute. 
One finds for the complex energy shift at first order: 

1 



7(1) 

'h,n£ 



h,nt 



2/i 



flu 



Vcc,h + Vcn,h 9 On Vnc,h 



1 47T 

27tT 



(2a 



1/2 3/2x 



+ 1 



nO / 1/2 







,3/2x2 



^no(O) 



6 



(4.3) 



where p* nQ is the cm. momentum of the neutral mesons after the decay of the bound state 
with quantum numbers (n,£ = 0). Taking into account that 



2 fj?a 3 
irn 3 



and introducing the isospin even and odd scattering lengths [Eqs. (3.9)], 



-(al /2 + 2ar), 



3/2>, 



1/ 1/2 3/2n 

3« ~ a o ), 



(4.4) 



(4.5) 



we obtain: 



E 



(i) 

h,n0 



h,nO 



-2/i' 



a 



n 



3 v a + a oJ> 



a" 



n 



3 v fl 



(4.6) 
(4.7) 
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The above formulas correspond to the expressions found by Deser et al. [1, 2, 3, 4] and 
provide the leading effects in the shift in the real part of the energy of the bound state 
and in the width of the decay into the neutral mesons. 

Numerical predictions about the scattering lengths are made from chiral perturbation 
theory (ChPT). They are summarized in Table 1. 



ChPT 


1/2 


3/2 


Tree [39] 
One-loop [40] 
Two-loop [41] 


0.14 
0.19 ±0.02 
0.220 


-0.07 
-0.05 ±0.02 
-0.047 



Table 1: Theoretical predictions for the 5"- wave scattering lengths from ChPT. 

On experimental grounds, the values of the scattering lengths are obtained from an 
extrapolation of high energy data down to the threshold. A complete analysis of the 
problem, using Roy [42] and Steiner [43] equations, has given [44]: 

m n al /2 = 0.224 ± 0.022, m n a /2 = -(0.448 ± 0.077) x 10" 1 . (4.8) 

Although the two-loop results of ChPT are not yet well understood [45, 46], one no- 
tices, with the improvement of the accuracy of calculations, a reasonable convergence of 
the theoretical estimates towards the experimental values. We shall adopt the central 
values of the experimental results (4.8) for the estimates that we shall do in the subse- 
quent sections for the various corrections to the lowest-order results; uncertainties in these 
corrections related to the central values will be neglected. For the present lowest-order 
results we obtain for the first three bound states the following estimates for the energies, 
decay widths and lifetimes (r), presented in Table 2. 



n i 


(eV) 


4 lj (eV) 


C (eV) 


7-W (10-^ S ) 


1 


-2898.61 


-8.86 


0.175 


3.76 


2 


-724.65 


-1.11 


0.022 


30.09 


2 1 


-724.65 


0.00 


0.000 


oo 



Table 2: Zeroth-order energies, first-order energy shifts, decay widths and lifetimes of 
hadronic origin of the first three bound states. 

In the following, we shall evaluate 0(a) corrections to these results. They originate 
from three effects: pure electromagnetic interaction, electromagnetic radiative corrections 
to the strong interaction and second-order perturbation theory effects of the bound state 
energy expansion. We shall consider them separately. 
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5 Electromagnetic interaction 



In this section we consider mainly corrections coming from pure electromagnetic inter- 
action. They arise in the channel tt~K + — > ir~K + from one- and two-photon exchange 
diagrams and also include vacuum polarization contribution. The corresponding dia- 
grams, together with the constraint diagram of the box-ladder diagram, are represented 
in Fig. 2. Vertex correction diagrams, associated with self-energy diagrams, contribute 
only at order a 5 In a -1 and may be ignored at the present level of precision (order a 4 ). 




x 



(b) 



(c) 



(d) 






(f) 



(g) 



00 



Figure 2: One- and two-photon exchange diagrams contributing to the pure electroma^ 
netic corrections. The diagram with a cross represents the constraint diagram. 



The on-mass shell one-photon exchange diagram (a) gives the contribution 

1 e 2 



cc, I7 



2,/s t 



2(s-mi 



m 



t 



(5-1] 



[e 2 = Ana], from which one has to subtract the Coulomb potential (3.14). An expansion 
of the total energy, according to Eq. (3.12), should also be made. 

In the category of two-photon exchange diagrams, we isolate the box-ladder, crossed- 
ladder and constraint diagrams, (b), (c), (d). The box- ladder and crossed-ladder diagrams 
have separately infrared divergences, as well as spurious singularities at threshold. To 
avoid their appearance, one must consider the sum of the above three diagrams, within 
which several mutual cancellations occur. The mechanism of cancellation is best under- 
stood with the threshold expansion method [34]. The leading part of potential momenta 
contribution, of order a 2 In a -1 , coming from the box-ladder diagram, is cancelled by 
that of the constraint diagram. The next-to- leading term, of order a 4 , vanishes in four 
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dimensions. Ultra-soft momenta do not contribute on the mass shell. Soft momenta 
contributions of box-ladder and crossed-ladder diagrams, of order a 3 In a -1 , cancel each 
other. One remains with 0(a 4 ) terms, the sum of which also vanishes on the mass shell. 
Therefore, the sum of the three diagrams (b), (c) and (d) does not contribute at order 
a 4 . Details can be found in the Appendix. 

The two-photon exchange diagrams have generally ultraviolet divergences. The coun- 
terterm lagrangian, which is a four-meson contact interaction term, does not contribute 
at order a 4 . 

Diagrams (e) and (f) contribute at order a 4 . The sum of their contribution is [Eq. 
(A.18)]: 



1 



Vcc,2, - g^= r=. (5.2) 



Diagram (g) does not contribute at order a 4 . 

The total corrective contribution of one- and two-photon exchange diagrams (without 
vacuum polarization), V cc ^ minus the Coulomb potential (3.14) and V CCj 2 7 , taking also 
into account the kinematic energy correction factor coming from the relativistic wave 
equation operator (last term of Eq. (3.24)), is: 

/Vo ij \ aA ^ (ca\ 

This coincides with similar results obtained from the Bethe-Salpeter equation (in the 
Coulomb gauge) [47] and from the Breit equation [26]. 

Another electromagnetic contribution is represented by the vacuum polarization dia- 
gram (h). Generally, in positronium, this diagram contributes at order 0(a 5 ). However, 
due to the mass differences between the electron, entering in the internal loop, and mesons, 
the contribution of the vacuum polarization diagram becomes enhanced [16]. It is numer- 
ically situated for the energy shift between 0(a 3 ) and 0(a 4 ) and turns out to be the 
most important correction to the strong interaction effect. It can be evaluated analyti- 
cally [13, 20]. We have evaluated it numerically for the first three bound states, from the 
expression of the corresponding local potential [48]. The results are compatible with the 
analytic evaluations and are presented below in Table 3. 

Apart from the pure electromagnetic corrections, we shall also evaluate here the strong 
interaction contributions to the meson electromagnetic form factors. The reason is that 
the interference effects between strong interaction and electromagnetism in the scattering 
amplitudes are evaluated in the literature without the above form factors, which are 
considered as parts of the one-photon exchange diagram (see Fig. 3). 

Defining the hadronic part of the electromagnetic form factor of pseudoscalar mesons 
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Figure 3: Strong interaction contribution to the electromagnetic form factors of mesons. 



as [37] 

F v {t) = l + Ur 2 )t + 0{t 2 ), (5.4) 
o 

where (r 2 ) defines the mean square charge radius of the meson, the corresponding contri- 
bution through the one-photon exchange diagram is: 

V ccMf = 2^(r 2 ). (5.5) 

[hff: hadronic form factor.] The sum of the contributions of the n and K mesons to the 
energy shift is: 

E$ fM = \^ ^ (<r V + (r 2 ) K+ ) 5 eo . (5.6) 

The mean square radii of the tt and K mesons were calculated from the data in Ref. 
[49] in the framework of ChPT to two loops. The values found there, which we use for 
the numerical evaluations, are the following: 

(r 2 )^ = (0.452 ± 0.013) fm 2 , (r 2 ) K+ = (0.363 ± 0.072) fm 2 . (5.7) 

Numerical estimates of the various electromagnetic contributions to the first three bound 
states are summarized in Table 3. 



n i 


Ell } +2h (eV) 


JSffl* (eV) 


< } / (eV) 


ESL (eV) 


1 


-0.147 


-2.561 


0.051 


-2.657 


2 


-0.025 


-0.296 


0.006 


-0.315 


2 1 


-0.006 


-0.025 


0.000 


-0.031 



Table 3: Electromagnetic corrections (elm) to the energy shift, composed of contributions 
of one- and two-photon exchanges ((1 + 2)7), vacuum polarization (vpol) and hadronic 
form factors (hff). 



18 



6 Strong interaction in the presence of 
electromagnetism 



To complete the evaluation of corrections at first order of perturbation theory, we have 
now to consider the interference effects between strong interaction and electromagnetism, 
including isospin symmetry breaking. Since we are interested by 0(a) corrections to the 
strong interaction effects found in Sec. 4, it is sufficient to consider diagrams with one 
photon propagator. The analysis is best carried out starting from the general relationship 
of the potential V with the two-meson irreducuble kernel K [Eq. (2.8)], which was already 
used in the pure strong interaction case. 

The kernel K can be separated into three parts, pure hadronic, Kh, pure electro- 
magnetic, Kry, and a part with interference between both, K^: K = + + K^. 
Replacing K in Eq. (2.8) and then linearizing with respect to K 1 + and subtracting 
the pure hadronic potential and the pure electromagnetic kernel K^, one ends up with 
the expression of the interference potential: 

V hl ={l-K h (Go - go))'' (tf 7 + K h7 )(l - (G - g ) K h )~ l - K T (6.1) 

Retaining the first few terms, we have for V^ 7 an expansion as follows: 

I4 7 = ^ 7 +^(Go-^o)^ 7 +^ 7 (Go-^o)^+^(Go-^o)^ 7 (Go-^o)^+--- • (6.2) 

In the present approximation, is the one-photon exchange kernel; also contains ef- 
fects of isospin symmetry breaking; internal propagators, such as Go, should be considered 
with physical masses. 

Typical diagrams, in the charged-charged (cc) channel, are represented in Fig. 4. 




(a) (b) (c) (d) 

Figure 4: Electromagnetic radiative corrections to the strong interaction in the charged- 
charged channel; diagram (d) is the constraint diagram of diagram (c); symmetric dia- 
grams to those, as well as self-energy diagrams, are not drawn. 
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Diagrams (a), (b), (c) represent one-photon exchanges in the t-channel, w-channel and 
s-channel respectively; diagram (d) is the constraint diagram of diagram (c). Also four 
other diagrams symmetric to the above ones exist; futhermore, one must also include the 
contributions of the self-energy diagrams. Each of those diagrams has infrared divergences 
on the mass shell; however, mutual cancellations occur by grouping several diagrams. The 
sum of the t-channel diagrams and the contributions of the self-energy diagrams is finite 
at order a 4 . The constraint diagram cancels the infrared divergence and spurious singu- 
larities of the s-channel diagram generated by the potential momenta [34] . The remaining 
part of the s-channel diagram, which is still infrared divergent and generated by the hard 
momenta, is associated with the w-channel diagram, which has similar singularities; their 
sum is free of divergences and of spurious singularities. Diagrams where the photon is 
emitted from the vertex (not drawn in Fig. 4) are not infrared singular and contribute as 
0(a 4 ). The sum of all the above diagrams is infrared finite and of order a 4 . It defines the 
regularized real part of the strong interaction vertex in the presence of electromagnetism; 
its value at the bound state energy differs from its value at threshold by an 0(a 2 ) term 
and therefore can be replaced by the value of the regularized real part of the vertex at 
threshold. (See Appendix for details.) 

In the charged-neutral (cn) or neutral-charged (nc) channels, t-channel and w-channel 
type diagrams are absent. In that case, one has to associate the part of the s-channel 
diagram generated by the hard momenta with the self-energy contributions of the external 
charged particles and the same type of cancellations as above operate, leading to the same 
qualitative result. 

The analysis done above can be repeated with more complicated diagrams. One has 
always to group several diagrams of the same class to reach, with the aid of the threshold 
expansion method [34], mutual cancellations of infrared divergences and spurious sin- 
gularities. The final results that we obtain are very similar to those found in Sec. 4, 
Eqs. (4.1), with the only difference that the pure strong interaction scattering amplitude 
is now replaced by the the regularized strong interaction amplitude in the presence of 
electromagnetism and isospin symmetry breaking: 



V cc,h+h-y — Vcc,h+h-y — Tie ^cc^h+h-y 



) 

thr. 



V —V = V = V =VpT re9 - 

v cn,h+h-y r nc,h+h-y y cn y nc cn 



= TZeTZ 9 ' 

thr. 



(6.3) 

thr. 

The threshold values of the regularized real parts of the strong interaction scattering 
amplitudes in the presence of electromagnetism deviate from the strong interaction S- 
wave scattering lengths [Eqs. (4.2)] by small amounts that we designate by (A(aQ+a>Q ))^ 7 
and (Aoq )h 7 - One then has: 

V cc ,h+h 7 = -47r((ao + % ) + (A(aJ + % )) /l7 J, 
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cn 



AnV2(a G + (Aa ) 

h-y ) i 



(6.4) 



where aj" and a are the pure strong interaction isospin even and odd scattering lengths, 
respectively [Eqs. (4.5)]. For further use, we define the following relative amounts: 



The various diagrams entering in the calculation of the strong interaction scattering 
amplitude in the presence of electromagnetism have also ultraviolet divergences. They are 
eliminated by the low energy constants of the effective chiral lagrangian in the presence 
of electromagnetism [50, 51]. 

Turning back to the infrared problem, we emphasize the following point. The con- 
straint diagrams have eliminated, to order a 4 , the infrared divergences and threshold 
singularities of the scattering amplitudes and allowed the definition of regularized real 
parts of them. In the literature, the scattering amplitudes are generally calculated in the 
scattering region; then, infrared divergences are eliminated by factorizing the Coulomb 
phase [52], which does not contribute to the cross section, and by combining the process 
under consideration with real soft photon emission processes [51, 53, 54]. The remaining 
threshold singularities are then subtracted to define a regularized scattering amplitude 
at threshold. In the present formalism, the use of constraint diagrams for the definition 
of the potentials circumvents the latter procedures on the mass shell (to order a 4 ) , pro- 
viding directly the regularized result. It can be checked explicitly that the pieces that 
are cancelled by the constraint diagrams are the same quantities that are subtracted in 
procedures dealing with the scattering amplitudes in the scattering region. Therefore, the 
regularized real parts of the amplitudes that we have defined in Eqs. (6.3) are identical 
to those defined in the literature. 

More complicated diagrams than those of Fig. 4 involve an increasing number of loops. 
However, in ChPT, the increase in the number of loops decreases the order of magnitude 
of the corresponding correction at low energy; this is why the chiral perturbation theory is 
organized in terms of the number of loops [36, 37, 38]. Therefore, the diagrams of the type 
of Fig. 4 represent the most important contributions to the interference effects between 
strong interaction and electromagnetism, apart from isospin partner mass difference in- 
sertion terms at the tree level, and practical calculation of these effects have been limited 
to them [51, 53, 54]. Next-to-leading effects are represented by diagrams of the type of 
Fig. 5. The latter diagrams contain infrared logarithmic contributions which partially 
enhance their order of magnitude and will be taken into account in Sec. 8 in conjunction 
with second-order perturbation theory effects. 

The evaluation for the 7rK atom of contributions of diagrams of the type of Fig. 4 and 
of mass shift insertion effects have been done by two groups of authors [53, 54]. Although 




6, 








(6.5) 
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(a) (b) (c) (d) 

Figure 5: Two-loop diagram and its constraint diagrams. 

their results do not exactly coincide, the numerical values that emerge for the corrective 
terms are close to each other. We use here the numerical values given in Ref. [54]. One 
has: 

€m = 0{e 2 ) + 0{e 2 p 2 )+0{\m u -m d \p 2 ) 

= (l.2- (0.3 ±3.2) + 0.2) x 10" 2 = (1.1 ±3.2) x 1(T 2 , (6.6) 

8% M = 0(e 2 )+0(\m u -m d \) + 0(e 2 p 2 )+0(\m u -m d \p 2 ) 

= (0.8 ± 0.5 - (0.8 ± 0.7) + (0.7 ± 0.2)) x 10~ 2 = (1.2 ± 0.7) x 10" 2 , 

(6.7) 

where 0(e 2 ) and 0( \m u — m d \ ) represent effects of electromagnetism and isospin symmetry 
breaking mass insertions at the tree level, 0(e 2 p 2 ) and 0( \m u — m d \ p 2 ) effects of electro- 
magnetism and isospin symmetry breaking at one-loop level; the uncertainties arise mainly 
from the electromagnetic low energy constants of the effective chiral lagrangian, which 
are estimated by their order of magnitude. 

The above modifications of threshold values of the amplitudes modify the results 
(4.6) and (4.7), obtained in the pure strong interaction case, with the inclusion of the 
corresponding corrections: 

E^ hj>n0 = -2/i 2 2j (a+ + a ")) (1 + O = £« (1 + O 

= <i (l±(l.l±3.2)xl0- 2 ), (6.8) 

(i) 2 q3 2 (i) (i) (1) 

= r« (l + (2.4±1.4)xl0- 2 ). (6.9) 
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7 Cancellation of divergences of three-dimensional 
diagrams 



We shall show in this section that divergences of constraint diagrams do not appear in 
physical quantities. Up to order a 4 effects, the only divergences that are introduced by 
constraint diagrams are those that cancel the infrared divergences of the on-mass shell 
scattering amplitudes. Other divergences that might appear in the perturbation series of 
the energy shifts are actually mutually cancelled. 

To observe the latter property, we go back to the matrix notation and express directly 
the scattering amplitude T' 2 [Eq. (2.21)], that defines the perturbation expansion of the 
energy shifts [Eq. (2.24)], in terms of two-particle irreducible kernels. Let K\ be the 
two-particle irreducible kernel of the unperturbed theory, that defines potential V\\ we 
have: 

^^(l-CGo-^i)" 1 . (7-1) 

Potential V 2 is defined as V 2 = V — Vi, where V is the the potential corresponding to 
the total interaction. It is related to the total scattering amplitude T [Eq. (2.3)] and to 
its two-particle irreducible kernel K [Eq. (2.8)], which in turn can be separated into two 
parts, after isolating in it K\\ 

K = K 1 + K 2 . (7.2) 

For the Green function G[ that enters in the definition of T' 2 [Eq. (2.21)], one can use a 
decomposition similar to that of Eq. (3.19): 

G[=g + g V igo + G';, (7.3) 

where G" corresponds to the part created by multiparticle or multiphoton exchanges. 

Replacing then V\ and V in the expression of T 2 in terms of K\ and K, respectively, 
and iterating T' 2 and G' l7 one obtains for T' 2 the following expansion, in which we have 
kept up to two propagator terms: 

% = K 2 + K 2 G K 2 + K 2 G Q K 2 G K 2 + K 2 G K 1 G K 2 
+K 2 G'[K 2 + %{G Q - g )K 2 + K 2 (G - g Q )% 
+K 1 (G - g )K 2 G K 2 + K 2 G K 2 (G - g^ 
+K X (G - s )Ki(Go - g )K 2 + K 2 (G - g )Ki(G - g^K x 

- g )K 2 (G - g^K, + ■■■ . (7.4) 

The important point to be noticed is that the effective propagator g is completely 
absent in terms containing only the kernel K 2 . g is present in terms that contain K\ 
on their left or right boundaries. Since the higher-order electromagnetic corrections with 
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respect to the Coulomb potential appear at order a 4 only at first order of perturbation 
theory, one may simplify the analysis by assuming that K\ corresponds essentially to 
the pure elecromagnetic part, A' 7 , and K 2 to the pure hadronic part, Kh, as well as to 
the interference part of the hadronic and electromagnetic interactions, Kh T The terms 
K1G0K2 or K2G0K1, in which K2 is represented by a constant contact term, have infrared 
divergences, the dominant part of which is cancelled by the constraint propagator g Q ; 
the non-dominant divergences are mutually cancelled, by considering together self-energy 
diagrams and eventually other combinations of t-channel and u-channel diagrams (cf. 
Sec. 6). The term K2G0K1G0K2, which in lowest order is represented by diagram (a) 
of Fig. 5, is finite and actually its constraint diagrams (b), (c) and (d) of Fig. 5 have 
been cancelled by similar terms that have appeared through the perturbative expansion 
of the bound state energy. Similarly, the constraint diagrams which were present in 
the definition of the hadronic potentials Vh (Sec. 3) also disappeared. Those diagrams, 
considered individually, contain ultraviolet or infrared divergences (linear or logarithmic). 
The present cancellation mecanism shows that their effect is irrelevant. However, the way 
of organizing the perturbation expansion in terms of the potentials V, rather than in terms 
of the kernels K or the amplitudes T, has the advantage, through cancellation effects by 
constraint diagrams, of naturally arriving at quantities that are smooth functions of the 
energy in the vicinity of the two-particle threshold, being continued, within the present 
approximations, to threshold (cf. Eqs. (4.1) and (6.3)). Actually, it is only the finite parts 
of their contributions that are relevant for that operation. This is why, we shall continue 
formulating the perturbative expansion of the bound state energy in terms of potentials, 
keeping in mind that in mutually cancelling constraint diagrams the same convention 
should be used when removing divergences. 

In summary, no ultraviolet or infrared divergences globally occur from three-dimen- 
sional diagrams in physical quantities up to oreder a 4 . 



8 Second order of perturbation theory 

We evaluate, in this section, the contributions of second-order perturbation theory effects, 
which are represented by the terms containing the subtracted Coulomb Green function 
G' c in Eq. (3.24). Here, the potentials that have significant contributions are the hadronic 
part and the vacuum polarization part of V cc . 

We first consider the hadronic part V cc ,h of V cc , which was already studied in Sec. 4. 
Its contribution to the energy shift at second-order of perturbation theory is represented 
by the following sum of terms: 



^^hh,rd — Pill Vcc,h "7^- Vcc,h + Vcc,h Vcn,h 9o n Vnc,h + V crlth g Qn Kc,/i 7^" Vcc,h <Pn£- (8.1) 
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(A negligible contribution, involving two g 0n s, has been omitted.) Since at the approx- 
imation we are working V CC) h and V cn ^ are constants in momentum space [Eqs. (4.2)], 
the contribution of G' c factorizes with its integrations, and the wave functions become 
projected on their values at the origin in x-space: 



d 3 p d 3 p' 1 , 



G' c (eS, P ,p') |^(0)| 2 V cc , h + 2V m , h g 0n V r 



(2tt) 3 (2tt) 3 2/x 



nc,h 



Defining 



(— ) 
K 2p' 

c( 2 ) 



III' 



d 3 p d 3 p' 1 (0) 



(2tt) 3 (2tt) 3 2p 



- 1/ / ^'c \ 

°M,n^ — v cc,h \ 2^ /ra£' 

and taking into account results (4.2), (4.3), (4.6) and (4.7) one obtains: 



£ 



(2) 



hh,n£ 



,(2) 



a( 2 ) f^w 



(1) 
h,n0 



which in turn yields: 
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hh,n0 
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/i/i,ra0' 



,(2) 
hh,n0 
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(2) 

hh,n0' 



(8.2) 

(8.3) 
(8.4) 

(8.5) 
(8.6) 

J h,n0 1 h,n0 

The calculation therefore amounts to that of the double integral of G' c /(2p). The 
latter is composed of three contributions [Eq. (3.19)]. 

The first corresponds to zero-photon exchange and its integral is equal to that of g 0c , 
i.e., to V / -6L( S )/(4 7 r) [Eq. (A.5)]. 

The second term corresponds to one-photon exchange; its integral is ultraviolet diver- 
gent, but this divergence is cancelled by that of the three constraint diagrams of Fig. 5 
(see Appendix for details). The finite part of the latter in turn cancels a finite logarithmic 
piece of the four-dimensional diagram (a) of Fig. 5; therefore, the finite part of the sum of 
the four diagrams of Fig. 5 becomes a smooth function. The finite part of the one-photon 
exchange part of the integral of G' c / (2/i) simply isolates the dominant logarithmic part 
of the four- dimensional diagram. 

The third term corresponds to the multiphoton exchanges and is finite. It can be 
calculated in several ways: either using an integration by parts in the variable p and 
isolating first the pole term to be subtracted [55], or integrating first with respect to the 
momenta and isolating at the end the pole term. 

The result is, for the finite part, using dimensional regularization and writing the 
contributions of the above three terms in successive order: 



,G' 



V 2/x 



(jf) 



pa 1 pa j . n , 
Ann 2tt a' 



+ 



pa 
2n~ 



•7) 
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where tj; is the logarithmic derivative of the Gamma function. 
Numerically, one finds, using Eq. (4.2): 



agL = 0.009, *g> = coos. 



We next consider the interference term between the hadronic and vacuum polarization 
parts of V cc . The corresponding energy shift is: 



Ql 

Z^hvpoiM - 2 W V cc,v P oi — V CC:h + V crhh g 0n V nCjh f ni , (8.9) 



Defining 



one has: 



which gives, using Eqs. (4.3): 

c(2) F (2) _ V(2) 

hvpol,n£ hvpol,n£ hvpol,nt 

= 2 (^ K c ,„ po ^) n r^(^£o-M^ )^ . (8.10) 
V ^ ' ^no(0) V 1 ' 

^hvpol.nO — (ptiO Vcc,vpol 7T^" ) + , . , (8-11) 
V W ^no(0) 

F m r (2) 

hvpol,nO nr(2) 1 hvpol,nO Ox( 2 ) (Q 1 0\ 

^(1) ~~ l0 hvpol,nOi (i) ~~ ^°hvpol,nO- 

^h,nO h,nO 

The above correction is finite. Using the expression of the vacuum polarization poten- 
tial V CC)Vpo i [48] and replacing G' c by a sum of contributions of intermediate states (discrete 
and continuous), one finds: 

*£Lio = °- 30a - °- 002 ' Cmo = 0.28a = 0.002. (8.13) 

On comparing the orders of magnitude of the individual contributions of V cc ^ and 
V cc ,vpoi from results (8.8) and (8.13), one deduces that the ratio of the contributions of 
the latter to the former is of the order of 1/4; this implies that the contribution of the 
quadratic piece in V CCtVpo i is in the ratio of 1/16 with respect to the quadratic piece of V cc ,h 
and hence can be neglected. 



9 Summary of results 

We calculated the main corrections to the lowest-order formulas of the energy shift and 
decay width of the nK atom. At lowest order, the energy shift and the decay width are 
given by the formulas obtained by Deser et al. (Sec. 4). They are expressed in terms of 
the iS-wave scattering lengths of the strong interaction nK — > ttK scattering amplitudes, 
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taken in the isospin symmetry limit; they are designated by E^ n0 and r^ , where n is 
the principal quantum number and i (here equal to zero) the orbital quantum number. 
The main corrections, of order a, that arise are the following: 

1) Pure electromagnetic interaction effects beyond the Coulomb potential, designated 
by E^i mnl (Sec. 5) in the bound state with quantum numbers n and t. They contribute 
only to the real energy shift. 

2) Electromagnetic radiative corrections as well as isospin symmetry breaking correc- 
tions to the strong interaction scattering amplitudes, the relative amounts of which with 
respect to the lowest-order results are designated by 5^ h and 5^ h , according to the 
charged-charged (n~K + — > tt~K + ) and chargrd-neutral (ir~K + — > tt K°) channels (Sec. 
6). 

3) Corrections coming from second-order perturbation theory of the expansion of the 
bound state energies (Sec. 8). They involve strong interaction type correlations, repre- 

(2) 

sented by the relative amount 5 hh n0 with respect to the lowest-order results, and strong 
interaction-vacuum polarization type correlations, represented by the relative amount 

X(2) 
hvpol,nO' 

The real energy shift, AE, and the decay width, T, including the 0(a) corrections, 
take the following expressions: 

AE n0 = ~2^ 2 ^ (<4 + a ) ( 1 + fi&fo + $hh,nO + 2 4«U,nO ) + E eL,nOi I 9 - 1 ) 

AE nl = E$l >nl , (9.2) 

TnO = SPnO/^^- ( fl f ( 1 + 2 ^cn,/ l7 + 2 $hh,n0 + 2 $hvpol,n0 ) ) ( 9 - 3 ) 

where the S"-wave scattering lengths and Oq are those of the pure strong interaction 
theory taken in the isospin symmetry limit, \x the reduced mass of tt~ and K + and p* n0 the 
cm. momentum of the neutral mesons after the decay of the bound state with quantum 
numbers (n, £ = 0). 

The numerical values of the various corrective terms for the first three bound states 
are summarized in Table 4. 



n t 


r(D 
u cc,h~f 


u cn, hf 


r(2) 

°hh 


r(2) 
u hvpol 


EZ (eV) 


1 


0.011 ±0.032 


0.012 ±0.007 


0.009 


0.002 


-2.66 


2 


0.011 ±0.032 


0.012 ±0.007 


0.008 


0.002 


-0.32 


2 1 


0.00 


0.00 


0.000 


0.000 


-0.03 



Table 4: Numerical values of the correction terms to the lowest-order formulas for the 
first three bound states. 
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The uncertainties come mainly from the low energy constants of ChPT in the presence 
of electrmagnetism which are poorly known and are taken into account through their order 
of magnitude. A more accurate evaluation of them, for instance by means of specific 
models, would considerably reduce those uncertainties. 

The decay width of the ground state and the energy splitting between the 2P and 2S 
states take the following forms in the presence of the 0(a) corrections: 

T 10 = 8p* w fi 2 a 3 (% f (1 + 0.046 ± 0.014), (9.4) 
(E 21 - E 20 ) = ~//a 3 (4 + a ) (1 + 0.023 ± 0.032) + 0.29 (eV). (9.5) 

These formulas allow one to extract from the experimental results on the decay width 
and the energy splitting the values of the strong interaction scattering lengths ag and ag~. 
We observe that the uncertainties in the corrective terms are much smaller in the decay 
width than in the energy splitting. This means that the decay width measurement will 
give us a more precise value for the scattering length tig than the measurement of the 
energy splitting for the combination (a^ + ). On the other hand, the latter quantity 
is sensitive, through its dependence on the low energy constants L4 and L$, to the Zweig 
rule violating effects [38] and therefore its precise knowledge is of crucial importance for 
the understanding of the chiral symmetry breaking mechanism and of the dynamical role 
of the strange quark [56, 57]. A precise knowledge of ag , which mainly depends on the 
low energy constant L5 [40], allow us to have a better insight into the ratio Fx I 'F n of the 
kaon and pion weak decay constants [38]. 
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A Three- and four-dimensional integrals 

Integrals are calculated with dimensional regularization, with dimension d close to 4. ~fi 
is the mass scale of the d-dimensional theory. 

For the various integrals with two, three or four propagators, we use a notation sim- 
ilar to that of Brown and Feynman [58]. In an elastic two-particle scattering process, 
we designate by p\ and p 2 the incoming particle momenta, with masses mi and rri2, re- 
spectively, and by p' x and p 2 the outgoing particle momenta, with the total momentum 
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P — (Pi + P2) = ip'i + p' 2 ) an d the momentum transfer q = ij>\ — p[) = (p' 2 — P2)', the 
Mandelstam variables are: s = P 2 , t = q 2 , u = (p\ — p' 2 ) 2 ■ 

For the propagators, we use simplified notations, where k is the loop momentum 
variable: 

1 ili 



(1) (pi — k) 2 — m\ + ie } (2) (p 2 + k) 2 — m 2 + ie'' 
1 ili 



^-l') (-P'i ~ k) 2 -m\ + ie' (-2') (-p' 2 + k) 2 - m 2 2 + is' 
1 -i 1 -i 



(0) k 2 + ie (3) (q-k) 2 + ie' 

The definitions of the integrals are: 

J = f [dk] - 1 F = / [dife] - ,.\ . . , = / [dfc] 



(a.i; 



1)(2)(0)(3)' 7 L J (l)(2)(3)' J 1 J (l)(2)(0)' 

G(1) = J [dk] mkwy G(2) = J [dk] (mwy A =I [dk] TWr (A - 2) 

where [dk] = ~p A ~ d d d k/(27r) d . Vector and tensor generalizations of these integrals corre- 
spond to the cases where momenta fc M or k^k u appear in the numerator of the integrands. 

Crossed diagrams involve integrals where (2) is replaced by (—2'). The corresponding 
integrals are defined as: 

J(l, -2') = / [dk] -—, — \, w r , (A.3) 
v ' J J [ 1 (l)(-2')(0)(3)' V ; 

etc. 

Constraint diagrams involve three-dimensional integrals which result from the s-chan- 
nel four-dimensional integrals by the replacement of the two propagators of the incoming 
particles by the single effective propagator g [Eq. (2.5)]. The corresponding integrals 
are, in the cm. frame: 

J° = -TT^F f [dk]5(k ) , F c = -^= f[dk}5(k 



2^ J L ^ U7 (l)(0)(3)' 2^J~s J ^ u '(l)(3)' 

BP = / [dk] S(k ) ^ A^-^j [dk] S(k ) i y . (A.4) 

The factor l/(2y/~s) that has been incorporated in the definitions takes account of the fact 
that constraint diagrams contain one more scattering amplitude than the four-dimensional 
integrals and the latter is defined with the factor i/(2y/s) [Eq. (2.2)]; the factor i is now 
contained in the definition of the meson propagator; the minus sign is reminiscent of a 
similar sign in front of g in Eq. (2.3). 
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The three-dimensional integral of the effective propagator g [Eqs. (2.4)-(2.5)] is, in 
the vicinity of d = 4: 



i2y/sA 



c 



1 

47T 



4-d 



d d - l k 



1 



k 2 + Z£ 



)[i + <H(-*i)-2 + M*) 



(A.5) 



'2 'V TW V 

where — &q(s) is taken positive and ijj is the logarithmic derivative of the Gamma function. 



Analytic continuation to positive values of &o(s) i s done with the replacement 

The integral A is: 
2 



A 



16vr 2 



d-A 



mim 2 . 
2 + In (- — T 



m\ — m 2 
2s 



In (^|) + Q(s) 



m 



(A.6) 



with 



Q{s) 



+ 



y/4s6q(s) f / s-(m 1 -m 2 ) 2 +y / s-{m 1 +m 2 )' 



y' s— (mi— rri2) 2 — y 7 s-(mi+m2) 2 



— Z7T 



(mi + m 2 ) 2 < s, 



(mi — W2) 2 < s < (mi + 7TI2) 2 , 



y4s (a) ^ / y(" -s+yH" 



y 7 (mi+m 2 ) 2 -s-y / (■m 1 -m 2 ) 2 



s < (mi — m,2) 2 . 



(AT) 



The integral of G , entering in Eqs. (2.6), (6.1) and (7.1), is equal to —i2y/sA. 
The integral J is equal, on the mass shell, to: 



J 



1 



167r tj-b 2 (s)s 



d-A 



X 



arctan 



m\ 



mi 



[8)8 



+ arctan 



m% 



[8)8 



(A.8) 



Its threshold expansion is obtained by taking into account the facts that \t\ <C s 
(mi + m 2 ) 2 , \bl(s) \ < s: 



J 



1 r 2 



16tt 2 t 



d-4 



7T 



3 a 



mim 2 



O^lncT 1 ), (A.9) 



where the order of magnitudes are evaluated with the counting rules of the QED bound 
states. The latter expansion can also be obtained with the threshold expansion method of 



Ref. [34]. The singularity in 1/ y— 6q(s)s is produced by the potential momenta, while the 
term in 1/ (m^j) is produced by the soft momenta. Ultrasoft momenta do not contribute 
at leading orders, while the hard momenta contribute at 0(a 3 ). 
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The expression of J(l, —2') is obtained from that of J by the replacement of the 



variable s by u. Its threshold expansion is obtained by noticing that u ~ (nil 
t ~ bl(s)(m 1 + m 2 ) 2 /(mim 2 ): 



m 2j 



.7(1,-2') 



16tt 2 i 



V>(1) + In i 



+ 0(a 3 lnaT 1 ). (A.10) 



The potential momenta do not contribute at leading order to J(l, —2') and the leading 
term comes now from the soft momenta. 

The three-dimensional integral J c is equal to: 



J 



c 



167T 2 t 



d-4 



-0(1) + In i 



47T7I" 



7T 



(A.H) 



s s 



Taking now the sum of the integrals J, J(l, —2') and J c , we find that J c cancels 
the potential momenta contribution of J, while the soft momenta contributions of J and 
J(l, —2') cancel each other: 

J + J(l,-2') + J c = 0(a 3 lna~ 1 ). (A.12) 

A similar type of analysis also applies to the sum of the three diagrams (b), (c) and 
(d) of Fig. 2 arising in QED. Here, the couplings being vector-like, one first decomposes 
the various integrals into a tensor sum, involving integrals of the types J, F, H, G, etc. 
and their vector and tensor associates. Taking into account the a 2 factor coming from 
the couplings, the sum of all these contributions vanishes up to order a 4 . (Notice that 
the constraint diagrams should be calculated as three-dimensional integrals involving g 
and one or two on-mass shell one-photon exchange diagrams.) 

The above result does not remain true in an off-mass shell formalism. Here, the 
scattering amplitude is no longer gauge invariant and slight differences arise. The sum 
(A.12) yields now a 0(a) term and generally the sum of the three previous diagrams 
behaves as 0(a> 3 ) [28]. The latter term, which is spurious, is cancelled by a higher-order 
diagram. Nevertheless, in the Fried- Yennie gauge [30], the same results as in the on- mass 
shell formalism occur. 

Integrals F and H appear also in electromagnetic radiative corrections to the strong 
interaction. Integral H appears in diagram (c) of Fig. 4; diagram (b) corresponds to its 
crossed diagram involving H(l, —2'); diagram (d) involves H ; diagram (a) involves H in 
the t-channel. H has the following threshold expansion: 



H 



7i 



327T 2 

1 

32tt 2 s 



r {s)s 
2 



d-4 



0(1) + In 



772 



d 



i-« 1 > +h w : 



+0(a 4 lna- 1 ; 



0i 



m\ — • 



02 



0i { 0i 

,2 



m\ — m 2 



1 , / 1 - 

V n( /v 

2 



(A.13) 
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where the dominant singularity comes from the potential momenta and the next-to-leading 
terms from the hard momenta. (The soft and ultrasoft momenta do not contribute at 
leading orders.) 

The crossed integral to H is: 



H(h-2') 



32tt 2 s 



d-4 



tp(l) + In 



AirJI 2 



+ 2 



1 



2 r 1 i , 1 n 1 w 1 > 



O^aUna- 1 ), (A.14) 



where only hard momenta contribute. The integral entering in the t-channel vertex func- 
tion (diagram (a) of Fig. 4) is H(l, — 1'; mi, mi): 



H(l, -1'; ,m 1 ,m 1 ] 



1 



327r 2 m? 



777 

^(l)+ln(^yj +0(a 4 lna- 1 ). (A.15) 



d-4 



It can also be obtained from the result (A.14), by taking in it the limit m 2 — > mi. 
if is equal to the opposite of the contribution of the potential momenta in H: 



H 



c 



TV 



32tt 2 



[8)8 



d-4 



(A.16) 



If E(p 2 , m 2 ) is the meson electromagnetic self-energy, taken for the moment with scalar 
couplings, then, after a mass-shell renormalization, it is the quantity ij^-| p 2 =m 2 that 
multiplies in lowest order the strong interaction vertex. This yields —1/2 of the value of 
the t-channel form factor at t — [Eq. (A.15)]. The self-energy contributions of the two 
external mesons of that form factor then cancel the latter completely at leading order. 
Similarly, the sum of H, H(l, —2') and H c yields a finite 0(a 3 ) term. 

The above results should be completed by incorporating the vector coupling of the 
photon. The latter does not change the leading behavior of H c . Concerning the four- 
dimensional integrals, the only modifications are through the hard momenta contributing 
with finite 0(a 3 ) effects. (Loop momenta in the numerators improve the infrared behavior 
of the object.) Taking also into account the factor a coming from the photon couplings, 
the sum of all contributions reduces to a finite 0(a A ) term. 

Diagrams where the photon is emitted from the vertex are not infrared singular and 
give contributions of order a 4 . 

In the case of the charged-neutral channel (process ir~K + — > 7r°K°), the integrals of 
the t-channel and u-channel form factors are absent. In that case, the cancellations occur 
between the sum H + H c and the contributions of the self-energies of the two external 
charged mesons, taking into account the vector coupling of the photon. The result is 
again a finite 0(a A ) term. 
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The integral that contributes to the dominant part of the constraint diagram (d) of 
Fig. 5 is denoted I cc ; it is: 

T cc 1 7-2(4-^) r d d 1 k d d l k' 



1^ J (2rr) d ~ 1 (2n) d ~ 1 b 2 (s)-k 2 + ie b 2 (s) - k' 2 + ie -(k-k') 2 + ie 

The integral that contributes to the dominant part of diagram (b) of Fig. 5, denoted 
I c is calculated in the following way. One first calculates the four- dimensional integral on 
the mass-shell; the latter is then integrated three- dimensionnally, by extending eventually 
the domain of validity of the momentum transfer squared t. The four-dimensional integral 
on the mass-shell is nothing but H [Eq. (A. 13)]. The latter is independent of t. Hence the 
three-dimensional integration reduces to that oiig , given by Eq. (A. 5) times i. The result 
is therefore the product of the right-hand sides of Eqs. (A. 13) and (A. 5), times i; since the 



integral of g is of order J — %(s), it is sufficient to retain the dominant singularity of H. 
One finds the opposite value of I [Eq. (A. 17)]. Hence the sum of the three constraint 
diagrams of Fig. 5 is given by the opposite value of l cc \ the finite logarithmic part of it 
cancels a similar term in the diagram (a) of Fig. 5 [34] . As we emphasized in Sees. 7 and 
8, the infinite part of I cc is cancelled by a similar term present in second-order of the 
perturbation theory expansion of the bound state energy. The finite part that has been 
retained in the right-hand side of Eq. (8.7), second term, corresponds to the contribution 
of the factor — 6q(s) of the logarithm. 

The integrals G [Eq. (A. 2)] are infrared finite and receive contributions at leading 
order from potential and soft momenta with the ratio (2)/(— 1): 

G (a) = = + 0(a 3 lna" 1 ), a = 1,2. (A.18) 

32m a v— t 
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